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Abstract 

The Laplace operator is considered for waveguides perturbed by a 



> 

o 
o 

. periodic structure consisting of congruent obstacles spanning the 

O I waveguide. Neumann boundary conditions are imposed on the pe- 

' riodic structure, and either Neumann or Dirichlet conditions on the 

. guide walls. It is proven that there are at least N (resp. A^— 1) trapped 

Q^l modes in the Neumann case (resp. Dirichlet case) under fairly general 

hypotheses, including the special case where the obstacles consist of 
line segments placed parallel to the waveguide walls. 



1 Introduction 
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■ Let TV be a positive integer. Consider the region in R : 



n = {{x,y):xe (-00, oo), y e (0, 2N)} - uZ=iOm, (1) 

with 

Om = {{x,y), X e [-a, a], y e [2m - 1 - g{x), 2m + 1 + g{x)]}. (2) 

Here a > 0, and g is a continuous function with g{x) G [0, 1) and g{±a) = 0. 
Thus fl can be viewed as a waveguide with N congruent obstacles placed 
periodically along the cross section. In 0, Linton and Mclvor studied the 
trapped modes in such regions under the hypothesis that g{x) was not iden- 
tically zero. Assuming Neumann boundary conditions on the obstacles, and 
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either Dirichlet or Neumann boundary conditions on the guide walls, they 
proved the existence of at least trapped modes in the Neumann case, and 
at least — 1 trapped modes in the Dirichlet case. The existence of these 
trapped modes was indicated earlier in numerical studies by Utsumomiya 
and Eatock Taylor and Evans and Porter These studies were moti- 
vated by a variety of possible applications to wave propagation in fluid and 
vibrating membranes; the reader is referred to 0] for a thorough discussion 
of these. 

Although the methods of Linton-Mclvor apply for a wide variety on as- 
sumptions of the geometry of the structure, they do not apply to the impor- 
tant special case where the structure consists of A^ identical line segments 
placed parallel to the guide walls, ie. g = 0. 

Furthermore, the numerical results in [Q, |^ also fail to indicate any 
trapped modes in this setting. The main purpose of this note is prove the 
existence of at least A^ (resp. A^ — 1) trapped modes in this setting for the 
Neumann (resp. Dirichlet) case. The methods of this paper also apply to the 
more general periodic structures described by Eqs. |I], 0, and this work might 
also be of interest because the upper bounds proven here on the associated 
frequencies will sometimes be sharper than those found in Linton-Mclvor. 

In addition, we consider the trapped modes of the region 

n={{x,y):xe (-00, oo), y G (0, 2N)} - U^-}0^, (3) 

with 

Om = {ix,y), X e[-a,a], y = 2m}. (4) 

For such regions with Neumann boundary conditions both on the obsta- 
cles and the guide walls, we prove the existence of A^ — 1 embedded eigen- 
values. The methods of this paper seem to fail in this case when one has 
Dirichlet boundary conditions on the guide walls. 

We note that a number of papers have proven the existence of at least 
one trapped mode for the case of a single line segment is placed parallel to 
the guide walls: see [@]) 01) and the references found therein. 

To prove the announced results, we use the symmetry of the problem to 
decompose the ambient Hilbert space into a direct sum of A^ + 1 invariant 
subspaces, as in . However, we differ in our choice of test functions. In the 
direction transversal to the guide walls, Linton and Mclvor's test function is 
essentially sinusoidal. In this paper, we choose a test function which has a 
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jump discontinuity across the obstacles, and which arises naturally from the 
Hilbert space decomposition. 



2 Statement of results and proofs 

We define the Laplace operator as 



A 



Let C^{Q) be the smooth functions of bounded support on Q. We shall 
study the self-adjoint operators Ai, A2 on -^^(fi), where Ai has operator 
core 

{uec^in):—u = o}, 

and A2 has operator core 

{u e C^{Q) : ^Uo, = 0, j = 1, . . . ,m, u\y=o = u\y=2N = 0}. 

Of course, Ai is simply the Neumann Laplacian, while we shall refer to A2 
as the Dirichlet case. 



2.1 Neumann boundary conditions 

In this section we prove the following: 

Theorem 1 Suppose Q satisfies Eqs. [J, Then, the operator Ai has at 
least N linearly independant trapped modes. Among the trapped modes, there 
exist N whose associated frequencies, denoted /i^, satisfy 

,rmT 2 

/im < (^) , m= l,...,iV. 

We begin the proof of the theorem by recalling the decomposition, used in 
[^], of L'^{fl) into Ai-invariant subspaces. Suppose / G L'^[0, 2N]. We extend 
/ to L2(— 00, 00) using the rules 

fiy) = f{-y), f{2N-y) = f{2N + y). 
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Set f{y) = Ei^=o fniv), where 

N N~l 

Uy) = ^ E cl/(2/ + 2n), (5) 

n=—N 

with 

c^,n = cos(mn7r/A^), (6) 

and 7^ = 2/(1 + 5^0 + ^mN)- Here is the Kronecker delta function. It is 
shown in that f{y) = J2m=ofmiy), and in fact 



N, 



with the Sm the image of the mapping f^fm- The subspaces {Sm} are mutu- 
ally orthogonal and are invariant under the Laplacian. We label Ail^^ = Am- 
Then 



info-ess (^r 



Am ■ 

Recall the Rayleigh quotient is given by 

where (j) is in the quadratic form domain of A^- For Neumann boundary 
conditions, the quadratic form domain of Am is 

Sm n H\n) = {ueSm- |Vm| g L\n)}. 

To prove the existence of eigenvalues below the essential spectrum of Am, 
(and hence the existence of trapped modes for Ai), it suffices to construct 
such that g(0) < 

Fix X and m, and for notational simplicity set = c„. We label 
the intervals (0, 1 - g{x)), (1 + g{x), 3 - g{x)), {2N - 1 + g{x), 2N) as 
/q, /i, . . . , /tv respectively. 

Let v{y) = 1 on Jq and elsewhere. Let v be the image of v under the 
mapping f^fm- Then, by Eq. |, 

N 

v{y) = ^cos{—-j)foTyeI,. (8) 
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Let b G [0, a), let a > 0, and define a piecewise different iable functions x 
and ipa on fl by 

0, |x| > a, 

1, |x| < b, 
xG(6,a) 

I a:G(-a,-6); 
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1, |x| < a. 
Our test function for Eq. ^ will be: 

,rmry. 



{x, y) = xix)v{y) + X^paix) cos( 



2N 



where A > is a parameter to be chosen later. Since v is in the image of 
the mapping f^fm, it follows that x'>^ ^ S^- Also, ipa{x) cos(^|^) G Sm by 
[^], Eqs.2.12, 2.13. Thus is in the quadratic form domain of Am- This test 
function can be compared to the one used in Eq.4.13. 
In what follows, it is convenient to set 

||v||^ = \vo\^ + 2\v2\^ + 2\v3\^ + ... + 2\vn-i\^ + \vn\^, 

where u,- = v\r.. 

Proposition 1 Let p = Then: 

>^'c.N + (II vi(l - g)iixr - P'X') - A^S^X sin(p(l - ^))) c^o; 
^ ^ A^iV/a + (A2iV(l - gix)) + ^ sin(if (1 -g))x+ || v||ix^(l - ^^(a:)) 1 dx 

here C = 4 /or m = 1, . . . , — 1 and C = 8 for m = N. 

The proof the this result appears in the appendix. 

We now complete the proof of the theorem. By Eq. ^ and the remarks 
that follow it, it suffices for the right hand side of the last equation to be 
negative. Note that the denominator is positive, and the same holds for the 
term J\x\<:aX{x) sin(p(l — g{x)))dx. Choose A sufficiently large that 

2^2 



(1 - g){Uf - p\') - A ^ / X sin(p(l - g)) < 0. 

x\<a mTC J\x\<a 
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Fixing this A, we then choose a > so that 

Ml / (1 - 9)iixT - pV) - A^-^ / X sin(p(l - g)) + X'aN < 0. 

J\x\<a mrC J\x\<a 

The theorem is proven. 

2.2 Dirichlet case 

In this section we prove the following: 

Theorem 2 Suppose Q satisfies Eqs. [J, Then, the operator A2 has at 
least N — 1 linearly independant trapped modes. Among the trapped modes, 
there exist N — 1 whose associated frequencies /ij satisfy 

^^m<{Y^)^ m = l,...,iV-l. 

For Dirichlet boundary conditions, we extend / G L^(0,2A^) to 00, cxd) 
via the equations 

f{-y) = -f{y), fm + y) = -f{2N - y). 

Then, using Eq. ^ as in the Neumann case, we have the decomposition 
L^{Q) = Q)m=oSm, where the Sm are mutually orthogonal and Sm are in- 
variant under A2. Setting Am = A2|s„, we have 



inf aess{Ar, 



TT^, m = 

:ff)2, m = l,...,iV. 



As in the Neumann case, we construct a test function for the Rayleigh quo- 
tient. In this case the quadratic form domain for Am is the closure in H^{Q)- 
norm of the set 

{ue Sm^ C^ifl) : support (m) n {y = 0} = support (m) n {y = 2N} = 0}. 
Let {Ij} be as in the Neumann case. Let 



1, ye Ij, 

0, elsewhere. 
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Denote by Vj the image of Vj under the mapping f^fm- Then, for j = 
1, . . . , iV — 1, we have by Eq. 

Vj{y) = ^(^i-^ ~ c-j~s), for yels 
mvr mir 

Here we have used cj = c_j, cj_2N = cj, and the identity cos(y4 — 5) — 
cos(A + B) = 2 sin A sin B. 

We chose j so that sin(^ ■ j) ^ 0. Thus Vj = if and only if m = or 
m = N . 

The test function is defined as: 

= x{'Av.j{y) + A^„(x) sin(^^), 

where A > is a parameter to be chosen later, and x, V'a as in the Neumann 
case. Note that xvj G S'm, and also that X'^j vanishes at ?/ = 0, ?/ = 2A^ (see 
Eq. I). Hence, x^i is in the quadratic form domain of Am- Also, by (P, 
Eq.2.15), ^/'q(x) sin(^l^) is in the quadratic form domain of A^-, and hence 
is in the quadratic form domain of A^- 

The proof of the theorem now follows from a word to word repetition of 
the argument used in the Neumann case. 



2.3 Line segments placed on ?/ = 2, 4, . . . , — 2 

In this section, we prove 

Theorem 3 Suppose Q satisfies Eqs. 0. Then, the operator Ai has at 
least N — 1 linearly independant trapped modes. Among the trapped modes, 
there exist N — 1 whose associated frequencies fim satisfy 

^^m<{^)^ m = l,...,iV-l. 

First, we note that the region f2 under these hypotheses satisfies the condi- 
tions necessary for the decomposition 

L\n) = So®...SN. 
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with f^fm defined exactly as above (see P, p. 3). 

We label the intervals (0, 2), (2, 4), . . . , {2N - 2, 2N) as Ji, /a, . . . , In re- 
spectively. Let {; = 1 on Ji, and v = elsewhere. 

Denote by v the image of v under the mapping f^fm- Then we have by 
Eq. I, 

v{y) = ^(c-i-s + c-s), for ye Is 

TTITT 1 TflTl 1 

= ^ cos(— ■ -) cos(— ■ (- + s)), for y G 

Here we have used cos(y4 — B) + cos(y4 + B) = 2 cos A cos B. 
Note that f = if and only if m = A^. 
The test function is defined as: 

y) = x{x)v{y) + Xtpaix) cos(^^), 

where A > is a parameter to be chosen later, and V'a as in the Neumann 
case. The proof of the theorem now follows by mimicking the argument used 
in the Neumann case. 



3 Appendix 

In this section we prove Proposition 1. 

We use the notation of Section 2.1. In what follows, we denote by {|x| < 
a} the set {{x,y) : |x| < a}. We have 

/ = 2X' r r cos\^)e-^-^^-^Uy dx 

J{\x\>a} Jx=aJy=0 Zl\ 

= \^N/a. 

Next, we calculate: 

J{\x\<a} •'{\A<a} 2N J{\x\<a} J {\x\<a} 2A' 

The first of the integrals on the right hand side of the last equation we 
compute as follows: 
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ly=2n-l+g{x) 2l\ J2N-l+g{x) 2l\ I 



n=l 



i=0 

= N I {l-g{x))dx. (10) 

J x=—a 

Here we have used the identity Z^ilo^ cos(^^^^) = 0, which follows from 

v^Af-l „2ijrmT/N _ n 
2^j=0 ^ ~ ^■ 

For the third integral on the right hand side of Eq. 10, a similar compu- 
tation yields 

r ,rmiy. 2N . ,nnt ,^ ... ,2mTci. . , 

j{|a;|<a} 2N Jx=-a mn 2N ~[ 2N 

ra 2N'^ mn 
= / Xsin( — (1-^)) ifm = l,...,(A^-l), 

Jx=-a mvr 2i\ 

= / Xsin( — (1 - (ix ifm = 0,A^. 

Jx=-a mn 2N 

Also, a straightforward calculation yields that 

pa 

v'^X^ = l|v||2 / - g{x))dx. 



f"- I , / smi — Tr-^j /x-^ ,mZm^^ \ , 



l{\x\<a} 

Thus, 

/ 101^ = A^iV/a+ r f A2iV(l - + ^ sin(^(l - ^?)) + ||v||^x'(l - g{x))] dx. 
Jn Jx=-a \ mn 2JS j 

(11) 

Here, C = 8 for m = 0, and C = 4 otherwise. 

In what follows, it is convenient to set ^ = P- Then 

W(x')' + 2At;cos(py)Cx')- (12) 

The fourth integrand is identically zero because ip'^x' = 0- Next, we calculate 
the first three terms in Eq. |T2|. First, 

AV / sm\py)^l = AV / (1 - cos\py)) + 2AV H [^^ e'^^^^-'^) sm\py) dy dx 

Jii J {\x\<a} J x=a J y=0 
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= AVA^ / (1 - g{x))dx + X'p'N/a; 

J x=—a 

here we have used Eq. |10|. Similarly, 

A^ / cos\py){ij'^f = X'Na, 



and 



Thus 

(13) 



|V0r = AViV / {l-g{x))dx+Xy'N/a+X'Na+\\v\\i / (x') 

Q J\x\<a J\x\<a 



The proposition now follows from Eqs. 13, 11 
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